home *** CD-ROM | disk | FTP | other *** search
/ CD ROM Paradise Collection 4 / CD ROM Paradise Collection 4 1995 Nov.iso / program / nrpas13.zip / GAMMQ.DEM < prev    next >
Text File  |  1991-04-29  |  621b  |  30 lines
  1. PROGRAM d6r7(input,output,dfile);
  2. (* driver for routine GAMMQ *)
  3. VAR
  4.    a,val,x : real;
  5.    i,nval : integer;
  6.    txt : string[25];
  7.    dfile : text;
  8.  
  9. (*$I MODFILE.PAS *)
  10. (*$I GAMMLN.PAS *)
  11.  
  12. (*$I GSER.PAS *)
  13.  
  14. (*$I GCF.PAS *)
  15.  
  16. (*$I GAMMQ.PAS *)
  17.  
  18. BEGIN
  19.    glopen(dfile,'fncval.dat');
  20.    REPEAT readln(dfile,txt) UNTIL (txt = 'Incomplete Gamma Function');
  21.    readln(dfile,nval);
  22.    writeln(txt);
  23.    writeln('a':4,'x':11,'actual':14,'gammq(a,x)':14);
  24.    FOR i := 1 to nval DO BEGIN
  25.       readln(dfile,a,x,val);
  26.       writeln(a:6:2,x:12:6,(1.0-val):12:6,gammq(a,x):12:6)
  27.    END;
  28.    close(dfile)
  29. END.
  30.